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Abstract 

We prove that the Hamilton Jacobi equation for an arbitrary Hamiltonian H (locally Lip- 
schitz but not necessarily convex) and fractional diffusion of order one (critical) has classical 
C 1,a solutions. The proof is achieved using a new Holder estimate for solutions of advection 
diffusion equations of order one with bounded vector fields that are not necessarily divergence 
free. 

1 Introduction 

We consider the Hamilton-Jacobi equation with fractional diffusion 

u t + H(Vu) + (-A) s u = 0. (1.1) 

Where H is a locally Lipschitz function. If the given initial data u(x, 0) = uq(x) is smooth and 
s > 1/2, the smoothing effect of the fractional Laplacian term is stronger than the effect of the 
nonlinear term in small scales, and it is well known that the solution u will remain smooth for 
positive time (see [5] and [3]), i.e. the problem is well posed in the classical sense. 

The case s = 1/2 is the most delicate because the two terms in the equation are of order one 
and their contributions are balanced at every scale. In this article we will show that the solution 
is also C 1,Q smooth in this critical case as well. 

For s < 1/2 we can show that the equation develops singularities (discontinuities of the deriva- 
tive) as a corollary of a construction of Kiselev, Nazarov and Shterenberg in [13] for the fractional 
Burgers equation. 

In [8], it was shown that equation (|1.1[) has a unique viscosity solution in W 1 ' 00 (i.e. Lipschitz) 
for any value of s € (0, 1) if the initial data u(—, 0) is Lipschitz. The equation was also studied in 
[9] and [11]. In [7], the following one dimensional equation is studied: 

v t + d x f(v) + (-A) s v = 0. 

This equation corresponds to the one dimensional case of (|1.1[) with / = H and v = u x . They 
prove that the equation has a smooth solution if s > 1/2, but they explicitly leave the case s = 1/2 
open, even in this one dimensional scenario. 

The equation (|1.1|) arises in problems of optimal control of processes with cn-stable noise. It is 
a particular case of the more general first order Isaacs equation 

u t - lu := u t - supinf (c ij + b ij -Vu+ f + . ^ ~ ^ a ij (y) dy) = 0. (1.2) 
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where i and j are two indexes ranging in arbitrary sets (controls), c lJ is a family of constant 
scalars, b lJ is a bounded family of vectors and the kernels a 13 {y) satisfy 

a lJ (y) = a lJ {— y) (symmetry) 
A < a 13 (y) < A (uniform ellipticity) 

The integral in (|1.2[) makes sense in the principal value sense (thanks to the symmetry assump- 
tion on the kernels a lJ ) as long as u G C ,a for some a > 0. The main result in this paper is that 
the solution u is C 1,a in both time and space, which means that the equations and (|1.2[) are 
well posed in the classical sense. 

Note that for any Lipschitz function H, the equation can be recovered from (|1.2p using 
a fixed constant a and writing H as 

H(p) = supinf(c y + b tJ ■ p). 

i 3 

Such representation of H always exists for some bounded family b lJ and d 3 if H is Lipschitz. If H 
is only locally Lipschitz, equation (jl.ll) takes the form (|1.2[) only a posteriori once we know that 
the solution u is Lipschitz, so that Vw stays in a bounded domain where H can be considered a 
Lipschitz function. 

In this paper we study problem (jl.2l) . The result holds for the particular case in every 
situation where the solution u is known to be Lipschitz. 

The idea of the proof is the following. If u is a solution of (jl.ip . any directional derivative 
v = u e would satisfy the linearized equation 

v t +w Vv + {-A) 1/2 v = (1.3) 

for w = DH(Vu). 

For a critical advection diffusion equation like (| 1 . 3[) . it was proved recently by Caffarelli and 
Vasseur [3] that the solutions become Holder continuous if divw = and w G BMO. Their 
method uses variational techniques in the style De Giorgi. Another proof was given by Kiselev 
and Nazarov in [T2]. In this article we establish a new proof that does not require w to be 
divergence free but it requires w G L°°. Our methods are non variational, as opposed to the De 
Giorgi style methods used in [4]. This new Holder estimate for an equation like (|1.3[) (Theorem 
16. lj) is interesting by itself and provides a non variational counterpart to the result in [3] and [12] . 
This is the key idea of this paper. Once we have it, the differentiability of the solution to (11.11) 
follows as a consequence. 

Note that we have no information a priori about the vector field w except that it is bounded 
by the Lipschitz norm of H: \w\ < A. Moreover, since w is not divergence free in general, we 
cannot make sense of (|1.3|) in the distributional sense. The only information we have is that 

v t -A\Vv\ + (-A) 1/2 v < 
v t + A\Vv\ + (-A) 1/2 v > 

which can be made sense of in the (Crandall-Lions) viscosity sense. So our Holder estimates 
depend on these two inequalities only. 

The key diminish of oscillation lemma is developed in section [S] The heart of the argument 
is Lemma \5. II In section [6l the results of section [5] are used to obtain the C 1,a estimates for u. 
In section [7] we discuss the case when the initial data is not Lipschitz. In sections [2] [3] and |4] we 
review the viscosity solution framework for Hamilton- Jacobi equations with fractional diffusion. 
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In section [51 we point out that in the supercritical case s < 1/2, the solution to can develop 
singularities as an immediate consequence of a construction of Kiselev and Nazarov [13] • 

The focus of this paper is on the C 1 '" regularity and not on the existence and uniqueness of 
a Lipschitz viscosity solution (which was proved in [5] for (jl.ip ). Nevertheless we sketch most of 
the necessary proofs in order to make the paper more self contained. 

2 Preliminaries 

We start by recalling the definition of the fractional Laplacian. The operator (— A) s is defined 
quickly using the Fourier transform as (— A) s u(£) = |£| 2s u(£). A more useful classical formula for 
the fractional Laplacian is 

(-A)X*) < X) ~:il +y) dy. 

Jr" \y\ n+2s 

The proof of this formula, as well as the computation of the precise constant C„ s , can be found 
in the book of Landkof ffj| . 

If a function u is C 1 '" and bounded, the operator (— A) 1//2 u is well defined and C a . This follows 
from the identity (—A) 1 / 2 ?/ = ^RidiU and the classical C a estimates for the Riesz transform. 
For a general elliptic operator of order one, the same statement holds and is proved in the following 
proposition. The important point to make is that once we prove that or (|1.2j) have a C 1 '" 
solution, this solution is classical. 

Proposition 2.1. Given any bounded symmetric function a < A and u G C ,a , the integro- 
differential operator 

f u{x + u) - u(x) 

LaU = l \y\ n+1 a{y)dy 

is a C a function whose C a norm depends only on A, and the dimension n. 

Consequently, any nonlinear operator of the form lu = inf j sup^ c lJ + b lJ ■ Vu + L a ij u is also 
C a if u G C ,Q as long as the family of vectors b 1 ^ is bounded and the kernels a l i are uniformly 
bounded. 

Proof. We will write the operator L a u as a classical singular integral operator applied to the 
gradient Vu. We write u{x + y) — u(x) as the integral of Vu • y along the segment from x to x + y 
and replace in the integral formula for L a u: 

f f 1 Vufx + sy) ■ y , . , 
Lau(x)= / \ J> y a(y)dsdy 

Ji" jo \y\ 

= / Vu(x + z) ■ dz 

where 

k(z) = J a (^—^j z ds. 

Since k is bounded and odd, then L a u € C a if Vu E C a by the classical Holder estimates for 
singular integrals. □ 

The theory of viscosity solutions developed by Crandall and Lions is very suitable to study 
solutions of equations (|1.1[) or (|1.2p . Here is a standard definition of viscosity solution adapted to 
the specific case of equation (|1.2p . 
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Definition 2.2. An upper (lower) semicontinuous function u is said to be a viscosity subsolution 
(supersolution) of (( 1 . 2|) if every time a C 2 function ip touches u from above (below) at a point 
(x, t) meaning that for some r > 



ip(x, t) = u(x, t) 
<p(y,s) > u(y,s) 



for all (y, s) such that \x — y\ < r and t — r < s < t. 



(resp. < u(y, s)) 



then if we construct the function v: 



v = 



( 



u 



in U :— B r (x) x (t — r, t + r) 
in (M™ x [0,+oo)) \ U 



then v satisfies 



Vt(x,t) 



lv(x,t) < 



(resp. > 0). 



A viscosity solution is a continuous function u that is at the same time a subsolution and a 
supersolution. 

This definition is the straight forward adaptation of the corresponding definitions in [T] and 
[5] for elliptic problems. 

The theory of viscosity solutions (See [5]) provides very general methods to obtain existence 
and uniqueness of solutions. Essentially, once a comparison principle is valid and barrier functions 
are constructed near the boundary of the equation, the existence follows by a standard Perron's 
method. 

Assuming that uo is continuous and bounded, cither system or (|1.2|) has a unique bounded 
continuous viscosity solution u. This was proved in [9] and [8] for (jl.ljl and Lipschitz initial data 
using a vanishing viscosity approximation. The same ideas could also be applied to (jl.2p . In the 
appendix we sketch a proof using Perron's method. 

3 Maximal Operators 

In this section we define the Pucci type extremal operators which are good over-estimators for the 
difference of two solutions of (|1.2[) . These operators were defined originally in [TB] and [3J- 

Definition 3.1. Given < A < A, we write the nonlocal maximal and minimal operators M + 
and M~ as 



where Su(x, y) := u(x + y) + u{x — y) — 2u(x) and x + and x stand for the positive and negative 
part of x respectively. 

These operators are the extremals of all uniformly elliptic integro differential operators of order 
one in the sense that for any function u which is C l a at the point x, 




M + u(x) 



sup 

a(y)=a{-y) 
^<a{y)<A 




a(y) dy. 



(3.1) 
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The equality above can be seen easily by averaging the value of y and — y in the integral. Since 

a{y) = a(-y), 

u{x + y)-u{x) If 5u(x,y) 

— W+ 1 — a{v) dy = 2 1 ~W Tr a[v) dy - 

Therefore, the choice of a which would make the integral larger is to choose a as large as possible 
where 5u(x, y) > and as small as possible where Su(x, y) < 0. 

A formula like (|3.1[) also holds for M~ replacing the sup by inf. Since M + and M~ arc 
a supremum and an infimum respectively of linear integro-differential operators with bounded 
kernels a, the result of Proposition 12.11 applies. Thus M + u and M~u will be C a if u £ C 1,a . 

We use the extremal operators to write an equation for the difference of two solutions of (|1.2[l . 

Lemma 3.2 (Equation for the difference of solutions). Let u be a viscosity subsolution and v a 
viscosity super solution of p.2[) . Assume \b lJ \ < A for all indexes i,j. then the function (u — v) is 
a viscosity subsolution of the equation 

{u - v) t - A\V(u -v)\- M + {u -v) <0. (3.2) 



On the other hand, if u be a viscosity supersolution and v a viscosity subsolution of (|1.2I) . 
(u — v) is a viscosity supersolution of 

(u - v) t + A\V(u -v)\- M~(u - v) > 0. 

In the case that the functions u and v are C 1,a , this lemma is straightforward from the definition 
of the operators M + and M~. When u and v are only semi continuous and satisfy the inequalities 
in the viscosity sense, the proof requires some work because the operators cannot be evaluated in 
the classical sense. These ideas are very well understood in the theory of viscosity solutions. We 
provide the proof in the appendix. 

Lemma 3.3 (Maximum principle). Let u be a bounded upper semicontinuous function which is a 
subsolution of the equation 

Ut - A\Vu\ - M + u < 0. 
Then for every t > 0, sup x u{x, t) < swp x u(x, 0). 

Proof. Let M > swp x t u(x,t). We do the proof by contradiction, assume that for some xq £ R" 
and to > 0, u(xq, to) > svcp x u(x, 0). Let H = M — u(xq, to) and e <C (u(xo, to) — svcp x u(x, t))/(l + 
A). 

Let g = 1/(1 + x 2 ) (any smooth bump function would suffice for this proof). Since g is smooth, 
M + g is continuous by proposition l2.il Moreover, we can rcscale it f(x) = g(Xx) with A sufficiently 
small so that M + f(x) = XM + g(Xx) < e and |V/(a;)| < e for all x £ M. n . Therefore the function 

<p(x, t) = M — hf(x - x ) + (1 + A)het 

is a supersolution of 

tp t - A\Vtp\ - M + ip > 0, 

for any value of h > 0. 

Note that for h = 0, ip > u and for any h > 0, <p(x, 0) > u(x, 0). Moreover, for h = 
H/(l — (1 + A)et), (p(xo,to) = u(x ,to) and for t > l/(e(l + A)), ip > u. Therefore there is a 
minimum value of h where ip does not stay strictly above u. For that value of h, <p will touch u 
from above at some point with t\ > 0. But this is impossible by the definition of viscosity 

solution since tp is a smooth supersolution of the equation. □ 
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Corollary 3.4 (Comparison principle). Let u be a viscosity subsolution and v a viscosity superso- 
lution of (|1.2[) . Assume u(x,0) < v(x,0) for every x S R™ and u and v are bounded. Then u < v 
in R" x [0,+oo). 

Proof. By Lcmma l3.2[ u — v satisfies the assumptions of Lemma 13.31 Thus u — v remains negative 
for all time and u < v. □ 



4 Lipschitz continuity 

If the initial data Uq is Lipschitz, we can provide a short proof that the same Lipschitz bound will 
be preserved by evolution of the equation (|1.2p using only the comparison principle. We show it 
in the next Lemma. 

Lemma 4.1. Assume u is a bounded viscosity solution of (|1.2|) and uq = u(— , 0) is a bounded 
Lipschitz function. Then u is also Lipschitz in x and for every t > 0, \\u(— , t)||j^_ < ll^olliip 

Proof. The Lipschitz constant ||uo||£i P = C is equivalent to the inequality 

u {x + y) < u (x) + C\y\ 

for every x and y. Now, for every fixed y, the function u(x + y, t) is a solution of (|1.2[) with initial 
data uq(x + y). On the other hand, the function u(x,t) + C\y\ is a solution of (|1.2[) with initial 
data uq(x) + C\y\. By comparison principle (Corollary 13. 4p . u(x + y,t) < u(x,t) + C\y\ for all 
x and all t. Since this argument can be repeated for all y, we obtain that u(—,t) is Lipschitz 
continuous, with Lipschitz constant C for alH > 0. □ 

Remark 4.2. The Lemma is based on comparison principle only, so the same proof can be applied 
to prove that viscosity solutions to remain Lipschitz continuous using only the comparison 
principle for 



5 The diminish of oscillation lemma 

In this section we prove the oscillation lemmas which will be used in section [6] to obtain Holder 
estimates. A somewhat simplified version of these Lemmas was used in [5]. 

The following is the key lemma of the paper, which provides a pointwise estimate from an 
estimate in measure. 

Lemma 5.1 (point estimate). Let u be an upper semi continuous function, u < 1 in W l x [—2, 0] 
which satisfies the following inequality in the viscosity sense in B2+2A x [— 2,0]. 

u t - A\Vu\ - M + u < e 

Assume also that 

\{u<0}n(B 1 x [-2,-l])| >(i. 

Then, if eo is small enough there is a 9 > such that u < 1 — 9 in B\ x [— 1, 0]. 
(the maximal value of £0 as well as the value of 9 depend only on A, A, A and n) 

Proof. We consider the following ODE and its solution m : [—2, 0] — > R: 



m(-2) = 0, 
m '(t) = c o\{x G B x : u(x,t) < 0}\ - Cim(i). 



(5.1) 
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The above ODE can be solved explicitly by the formula 



m(t) = J c \{x : u(x,s) < 0} n B x | e ~ Cl(t ~ s) ds. 

The strategy of the proof is to show that if Co is small and C\ is large, then u < 1 — m{t) + 2eo 
in B x x [-1,0]. Since for t G [-1,0], 

m(t) > c e- 2Cl \{u < 0} nBi x [-2,-l]| > c e~ 2Cl \i. 

We can set 6 = c$e~ 2Cl fx/2 for £o small and obtain the result of the Lemma. 

Let /3 : R — > R be a fixed smooth non increasing function such that /3(x) = 1 if x < 1 and 
/3(x) = if x > 2. 

Let 6(x, i) = /3(|x| + At) = (3(\x\ — A\t\). As a function of x, b(x, t) looks like a bump function 
for every fixed t. At those points where 6 = (precisely where \x\ > 2 — At = 2 + A\t\), M~b > 0. 
Since b is smooth, M~b is continuous by Proposition ^ . 1 1 and it remains positive for b small enough. 
Thus, there is some constant f3\ such that M~b > if b < 

Assume that it(x, i) > 1 — m(t) + Eq(2 + t) for some point (x, t) E B\ X [—1, 0]. We will arrive 
to a contradiction by looking at the maximum of the function 

w(x, t) — u(x, t) + m(t)b(x, t) — £o(2 + t). 

We are assuming that there is one point in B\ x [—1,0] where w(x,t) > 1. Let (xo,to) be the 
point that realizes the maximum of w: 

w(xo,to) = max w(x,t). 

R»x[-2,0] 

Note that this maximum is larger than one, and thus it must be realized in the support of b. So 
|xo| < 2 + A|to| < 2 + 2A, and (xo, to) belongs to the domain of the equation. 

Let <p(x, t) := w(xo, to) — m(t)b(x, t) + £o(2 + 1). Since w realizes its maximum at (xo, to) then 
if touches u from above at the point (xo,to)- We can then use the definition of viscosity solution. 
For any neighborhood U of xq wc define. 



v(x, t) 

and wc have 



u(x, t) if x $l U 
ip(x, t) if x G U. 



v t - A\Vv\ - M + v < eo at (x ,t ). (5.2) 
We start by estimating i> t (xo,£o)- 



v t (x Q ,to) = -m'(t )b(xo,t ) - m(t )dtb{x Q ,to) + e 

= -m'(t )b(x ,t ) + m(to)A\Vb(x ,t )\ + e 



(5.3) 



We also note that 

|V«(« ,to)| = m(<)|V6(x ,t )|. (5.4) 

The delicate part of the argument is to estimate M + v correctly. Let us choose U to be a tiny 
ball B r for some r<Cl, The computations below are for fixed t = to, so we omit writing the time 
t in order the keep the computations cleaner. 

Since u + mb attains its maximum at xo, <5w(xo, y) < — m<5&(xo,y) for all y. Replacing this 
inequality in the formula for M + and M~ we can easily obtain M + v(xq) < — m M~b{xo), but 
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this is not sharp enough. We need our estimate to take into account the measure of the set 

{u<o}nBi. 

Let y € M. n be such that u(xo + y) < 0. We estimate Su(xo, y) + m 5b(xo,y). 

Su(x ,y) + m Sb(x Ql y) = 

u{x + y) + m b(x Q + y) + u(x - y) + m b(x Q - y) - 2u(x ) - 2m b(x ) 

Since u + mb attains its maximum at xq, 

< u(x + y) + m b(x + y) — u(x n ) - m b(x ) 
Since u(x n ) + m b(x ) = w(xo 7 t ) + £o(l + to) > 1, 

< + m - 1 
We choose Co small so that m < 1/2 and 

Su(x + y) + m Sb(x + y) < -- 

Now we estimate M + v(xo,to), we start writing the integral 

+ m(t ) f ASb{xo,y)- - XSb(x ,y) + f ASu(x , y) + - X5u(x , y) _ 

M^v{x ,to) = — / — — dy+ / — — dy 

2 J Br \y\ + JM. n \B r \y\ + 

We estimate 8u(xq, y) by above by —Sb(xo,y) except at those points where xq + y is in the good 
set G :— {u < 0} Pi £>i where we use that Su + mSb < —1/2 

. ur ,, + \,f A5(u + mb){x Q ,y) + - XS(u + mb)(x , y)~ 
<-m{t )M b{x ,t )+ / — — dy 

JG\s r w\ 

< -m{t )M-b(x ,t ) + f dy 

J g\b t \y\ 

< -m(t )M-b(x ,to)-c \G\B r \ 

for some universal constant cq (this is how Co is chosen in (|5.1[1 ). Note that as r — > 0, the measure 
of the set \G \ B r \ becomes arbitrarily close to \G\ = {x 6 B x : u(x,t ) < 0}. 

We consider two cases and obtain a contradiction in both. Either b(xQ,t Q ) > j3\ or b(xo,to) < 

Pi. 

Let us start with the latter. If b(xo,to) < /3i, then M~b(xo,to) > 0, then 

M+v{x ,t a ) <-c \G\B r \ (5.5) 
Replacing (jO]) . and (|53|) into ([5~2l we obtain 

e > -m'(to)&(ar ,to) + + c |G \ B r \ 

but for any C\ > this will be a contradiction with (|5.ip by taking r small enough. 

Let us now analyze the case b(xo, to) > fi\. Since b is a smooth, compactly supported function, 
there is some constant C (depending on A), such that |M~&| < C. Then we have the bound 

M + v{x ,t ) < -m(t )M-b(x ,t ) + c \G\B r \ < -Cm(t ) + c \G \ B r \ 
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Therefore, replacing in (|5.2[) . we obtain 

£q > -m'{t )b(x ,t ) +e + c \G\B r \ - Cm(t ) 

and we have 

-m'(t )b(x ,t ) -Cm(t ) + co\G\B r \ < 0. 
We replace the value of m'(to) in the above inequality using (|5.1[) and let r — s> to obtain 

(Ci6(x 0l to) - C)m{t ) + cq(1 - 6(0:0, to))|G| < 0. 
Recalling that 6(xo,to) > we arrive at a contradiction if Ci is chosen large enough. □ 
In the next lemma we use the notation Q r to denote the cylinder 

Q r := B r x [— r, 0]. 

We prefer to use this notation instead of C r because we reserve the letter C for constants. 

Lemma 5.2 (diminish of oscillation). Let u be a bounded continuous function which satisfies the 
following two inequalities in the viscosity sense in Qi 

u t - A\Vu\ - M + u < (5.6) 
u t + A\Vu\ - M~u > (5.7) 

There are universal constants 9 > and a > (depending only on A, the dimension n, and the 
ellipticity constants X and A ) such that if 

\u\<l in Qx := B\ X [—1, 0] 

\u(x,t)\ < 2\(4 + AA)x\ a - 1 in (R n \Bi) x [-1,0] 

then 

osc u < (1 — ff) 

Ql/(4+4A) 

Proof. We consider the rescaled version of u: 

u(x, t) = u((4 + AA)x, (4 + 4A)t). 

The function u will stay either positive of negative in half of the points in B\ x [—2,-1] (in 
measure). More precisely, either {u < 0}n(Bi x [-2, -1]) > \Bi\/2 or {u > 0}n(Bi x [-2, -1]) > 
|Bi|/2. Let us assume the former, otherwise we can repeat the proof with —it instead of u. 

We will conclude the proof as soon as we can apply Lemma 15.11 to u. The only hypothesis 
we are missing is that u is not bounded above by 1. So we have to consider v = min(l,u) and 
estimate the error in the right hand side of the equation. We prove that if a is small enough, then 
v satisfies 

v t - A\Vv\ - M + v < e 

for a small e and we can apply Lemma |5. II 

Note that inside Q4+4A, u < 1, thus v = u. The error in the equation in Q2+2A comes only 
from the tails of the integrals in the computation of M + v. Indeed, if ip touches v from above at 
a point (x,t) £ then it also touches u at the same point. Choosing a small neighborhood 

U of x and constructing 
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I w in U \ Lp in U 

tt'l = < W2 = \ 

I u outside U \v outside U 

We see that dtwi — A|Vifi| — M + w\ < at (x,t). On the other hand dtW\ — A\Vw\\ 
dtW2 — AlVu^l- So we estimate M + wi — M + ui2- 

M + Wl (x,t)-M + w 2 (x,t)= [ A(u(x + y)-l) + j: j^ T 



< 



x+y^Bi+iA \U\ 

A(2((4 + 4A) 2 |x + 2/ |r-2)+^ T < £0 

4+4A \y\ 



if a is chosen small enough. 

Thus, for any test function (p, dtW2 — A|Vt«2| — M + W2 < £o and u satisfies that inequality in 
the viscosity sense. We can then apply Lemma \b. II to v and conclude the proof. □ 



6 C 1,a regularity 

Theorem 6.1 (Holder continuity for advection-diffusion equations). Let u be a bounded continuous 
function which satisfies the following two inequalities in the viscosity sense in l n x [0,t] 

u t - A\Vu\ - M+u < (6.1) 
u t + A\Vu\ - M~u > (6.2) 

There there is an a > ( depending only on A, the dimension n, and the ellipticity constants X 
and A ) such that for every t > the function u is C a . Moreover we have the estimate 

\x — v\ a + \t — s\ a 

\u(x,t) -u(y,s)\ < C||u(-,0)|| ioo 1 for every x,y el" and < s < t 

(6.3) 

Equivalently we can write the estimate as, 



C 



ll M llo(R"x[t/2,t]) < ll«(-.0)|| 



L ' 



Proof. For any (xo,io), we consider the normalized function 

v ( x > *) = Trnr — u ( x o + to%, h(t + l)). 

We prove the C a estimate (|6.3|) at every point (xq, t ) by proving a C a estimate for v at (0, 0). 
Note that since the L°° norm of u is non increasing in time, then \v\ < 1. Moreover, u is also a 
solution of the same equation (|6.1|) and f|6 . 2[) . 

Let r = 1/(4 + AA). The estimate follows as soon as we can prove 

osc v < 2r ak for k = 0, 1, 2, 3, . . . (6.4) 

We will prove (|6.4p by constructing two sequences afc and such that ak < u < frfc in <5 r fc, 
frfc ^ a k — 2r Qfc , afc is nondecreasing and 6^ is non increasing. We will construct the sequence 
inductively. 
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Since \v\ < 1 everywhere, we can start by choosing some ao < inf v and bo > supt; so that 
bo — a o = 2. Assume we have constructed the sequences up to some value of k and let us find afe+i 
and b k +i- 

We scale again by considering 



w(x, t) = (v(r k x, r k t) - {a k + b k )/2)r 



—ak 



Therefore we have 



\w\ < 1 in Qi 

\w\ < 2r~ ak - 1 in Q r - k therefore \w{x, t)\ < 2\r~ 1 x\ a - 1 where \x\ > 1. 



If a is small enough, we can apply Lemma 15.21 to obtain oscg r < 1 — 9. So, if a is chosen 
smaller than the a of Lemma (|5.2[) and also so that 1 — 8 < r a , we have oscQ r < r a , which means 
0SC Q k +i — r a ( k+is> so we can find a k +i and b k +i and we finish the proof. □ 



The following theorem is actually a corollary of Theorem 16.11 but we state it as a theorem 
since it is the main result of this paper. 

Theorem 6.2 (C ,Q regularity for HJ with critical diffusion). Let u be a bounded continuous 
function which solves (|1.2[) in R™ x [0,i\. Assume the initial data u(— ,0) is Lipschitz continuous. 
There is an a > ( depending only on A, the dimension n, and the ellipticity constants A and A ) 
such that for every t > the function u is C > a in both x and t. Moreover 

C 

IMIc^(R"x[tAt]) < ^ IK- 0)IIl; p 



Proof. The proof follows by applying Theorem 16.11 to incremental quotients of u. We start by 
proving the regularity in the space variable x. 
For any vector e, the incremental quotient 

. . u(x + e,t) — u(x, t) 

V e {x,t) = p- 

|e| 

is bounded in L°° by ||u|| iip and satisfies the hypothesis of Theorem 16. II Therefore 

C 

Kllo(R»x[t/2,t]) < ^ IK->0)|| L . tp . 

uniformly in e. Thus V^w is C a with an estimate 

C 
f 



Since u is C ,a in space, the operator right hand side in the equation (jl.2[) . Iu is bounded (and 
Holder continuous). Therefore u t is bounded. Therefore, we can consider an incremental quotient 
in time 

. . u(x, t + h) — u(x, t) 
Wh(X,t) = , 

and Wh will be bounded in L°° independently of h. Moreover, uih satisfies the hypothesis of 
Theorem 16. 1[ then Wh is C a independently of h, which implies that ut is Holder continuous as 
well with an estimate 

C 

IMIc»(R"x[t/2,t]) < ^ ll«(-.0)lliip. 

which finishes the proof. □ 
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7 Non smooth initial data 



We have proved in Theorem 16 . 21 that if the initial data u(— , 0) = uq is Lipschitz, then the solutions 
u immediately becomes C 1,a for t > 0. In this section we will show that the Lipschitz condition 
on uq is not necessary. 

For equation (|l.ip the situation is somewhat different. If H is globally Lipschitz, then it is 
a particular case of (|1.2[) and the following theorem applies. If H is only locally Lipschitz, then 
we must first show that u is Lipschitz in order to apply our theorems for (|1.2[) . and therefore the 
following result is not relevant. 

Theorem 7.1 (C ' a regularity for non Lipschitz initial data). Let u be a bounded continuous 
function which solves (|1.2[) in R" x [0,t]. There is an a > (depending only on A, the dimension 
n, and the ellipticity constants X and A) such that for every t > the function u is C ' a in both 
x and t. Moreover 

C 

\\ u \\c 1 -"{R"x[t/2,t]) - 0)IIl= 

The proof of Theorem 17.11 uses the following Lemma from [2] (section 5.3) in order to improve 
the regularity estimates on u applying Theorem 16.11 repeatedly. 

Lemma 7.2. Let < a < 1, < /3 < 1 and K > be constants. Let u E L°°(R n ) satisfy 
IM|_loo — Define, for h E W 1 , with < \h\ < 1, the incremental quotient 

u(x + h) - u(x) _ 
vp, h (x) = — r^p — for x £ R . 

Assume that vp t h £ C a (W n ) with ||u^,ft|| CQ ( R „j < K for any < \h\ < 1. We then have 

1. If a + 13 < 1 then u E C Q+/3 (M") and \\u\\ ca+p < CK. 

2. If a + [3 > 1 then u <E Lip(R" ) and \\u\\ Lip < CK. 
where the constant C depends only on a and j3. 

We use the Lemma above in order to probe Theorem 17.11 



Proof of Theorem \7.1\ All we must prove is that the function u{— ,t) becomes Lipschitz for any 
t > 0. Then we can apply Theorem 16.21 to obtain the C 1,a estimate. 

Let Co = 10 (the operator I applied to the constant zero function, which returns a constant). 
We have that v(x,t) = cat is a particular solution to ()1.2|) . We apply Theorem 16.11 to u — v to 
obtain that u — v is C a for t > 0. Therefore u becomes C a for t > 0. 

Now, we apply iteratively Theorem 16. II to incremental quotients of u of the form 



. . u(x + h) — u(x) 
V0,h{*) = ■ 

We start with f3 = a. Since u G C a , then vp.h is bounded in L°° for any t > independently of 
h. From Theorem 16. 1[ we have that Vp^ becomes uniformly C a for t > 0. But then from Lemma 
17.21 u E C 2a . We repeat this procedure to obtain u E C 3 for f3 = 2a, 3a, 4a, . . . , until we reach 
the fcth step when kf3 > 1 and we obtain u E Thus we can apply Theorem l6.2l and finish the 
proof. □ 

Remark 7.3. Under some special assumptions on H, it may be possible to prove that solutions 
of p.ip whose initial data is only uniformly continuous become Lipschitz for t > 0. This is well 
known for the classical Hamilton- Jacobi equation without fractional diffusion (see |15|). 
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8 The supercritical case. Non differentiability. 



We can make an example of a non differentiablc solution of (|1.2[) with smooth initial data for any 
s < 1/2 in one dimension for H(p) = \p\ 2 as a corollary of a result of Kiselev and Nazarov [13]. 

Theorem 8.1. For any s < 1/2 there is one smooth function uq : R — > R swc/i i/iai t/ie equation 

ut + ^ + (-A) s « = 

does not have a global in time C 1 solution with u{— ,0) = «o- 
Proof. The derivative w = Wj satisfies the fractional Burgers equation 

v t + w x + (-A) s w = 0. 

In [13], it is shown that this equation can develop shocks for any s < 1/2. □ 

Appendix 

Proof of Lemma 13.21 

In this appendix we prove Lemma 13.21 We use Jensen's sup-convolutions idea jlO) , which is a 
standard method in the theory of viscosity solutions. 

The definition of viscosity solution satisfies the following important stability condition with 
respect to half relaxed limits (See [3] or [6]). 

Proposition 8.2. Let v k be a sequence of viscosity supersolutions of (jl.2p . Let v* be the lower 
half relaxed limit (equivalently the Y -limit): 

v»(x,t) = liminf v k {y, s). 

fe— > + oo 
3/->x 
s— >i 

Then v* is also a supersolution of (jl.2[) . 

If Uk be a sequence of viscosity subsolutions of (|1.2[) . Let u* be the upper half relaxed limit 
(equivalently —u* is the T -limit of —Uk): 

u*(x,t) = limsupu/ c (?/, s). 

k— S- + OG 

y— * x 

s— H 

Then u* is also a subsolution of (| 1 . 2[) . 

We start by recalling the definition of inf and sup convolutions. 
Definition 8.3. Given an upper semicontinuous function u, the sup-convolution u e is defined as 

u E (x,t)= sup u(y,s)--(\x-y\ 2 + \t-s\ 2 ). 

y£R",sG[0,+oo) £ 

Conversely, for a lower semicontinuous function v, the inf- convolution v e is defined as 
v e (x,t)= inf v(y,s) + -{\x-y\ 2 + \t-s\ 2 ). 

yGR™,sS[0,+oo) £ 
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The importance of this construction is that it regularizes the function u making it semiconvex 
in the case of u £ and semi concave in the case u £ , and therefore Lipschitz, while preserving the 
condition of sub or supersolution as it is pointed out by the following standard lemma. 

Lemma 8.4. Ifu(x,t) is a viscosity subsolution of (jl.2[) . then so is u £ . 
Ifv(x,t) is a viscosity supersolution of (|1.2|) . then so is v s . 

Proof. The proof follows by observing that u £ is a supremum of translations of u and the equation 
p.2[) is translation invariant. The same idea applies to v s . □ 

In the next proposition we point out, without a proof, the standard properties of inf and sup 
convolution (See for example [B] for the proofs and further discussion). 

Proposition 8.5. • The function u £ is semiconvex in the sense that every point has a tangent 
paraboloid from below with opening 1/e. 

• The function v e is semi concave in the sense that every point has a tangent paraboloid from 
above with opening 1/e. 

• The functions u e and v £ are Lipschitz, with Lipschitz constant bounded by | 

• u £ converges to u and v £ converges to v in the half relaxed sense (i.e. v £ T -converges to v 
and —u £ IT '-converges to —u) as e — > 0. 

Now we are ready to write the proof of Lemma 13.21 in the same fashion as the corresponding 
theorem in [3]. 

Proof of Lemma \3.2l By Lemma 18.41 u £ is also a subsolution of (jl.2j) and v e a supersolution. 
Moreover, by Proposition 18. 5[ u £ — > u and v s — > v in the half relaxed sense. By the stability of 
viscosity solutions under half relaxed limits, it is enough to show that w = (u £ — v £ ) is a subsolution 
of dO l. 

Let ip be a smooth function touching (u £ — v e ) from above at the point (x,t) (with t > 0). 
Note that for any e > 0, both functions u £ and —v e are semiconvex, which means that they have 
a tangent paraboloid from below of opening Since the smooth function ip touches (u £ — v £ ) 
by above at the point (x,t), then both u £ and — v E must be C 1 ' 1 at the point (x,t) (meaning that 
there are tangent paraboloids from both sides). But this means that dtu* , dtv*, Vit* and Vi>* are 
well defined classically at the point (x,t). Moreover it also means that we can evaluate lu £ (x,t) 
and lv E (x,t) in the classical sense. But then we compute directly at the point (x, t), 

(u £ - v £ ) t - A\\7(u £ - v £ ) \ - M+(u £ - v e ) < d t u* - lu* - d t v* + Iu, < 

which clearly implies that also 

ipt(x,t) - A\Vip(x,t)\ - M + tp(x,t) < 

since at the point (x,t), ipt = (u £ — v £ ) t , V<p = V(it £ — v £ ) and M + ip > M + (u £ — v £ ), all in the 
classical sense. 

This proves that u — v satisfies (|3.2|) . The other inequality follows in a similar way. □ 
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Comparison principle for (II. ip 

Given a locally Lipschitz kernel H, we can apply Lemma [5721 to a solution of only a posteriori 
if we know that the viscosity solution u is Lipschitz, so that we can restrict H to a bounded domain. 
The Lipschitz continuity of the solution u of (|1.1[) is a simple consequence of comparison principle 
if the initial data is already Lipschitz (Lemma 14.11) . However, we need to prove the comparison 
principle for equation (jl.ll) without using Lemma 13.21 This can be done by applying the idea of 
the proof of Lemma 13.31 directly to the difference of sup and inf convolutions. However one must 
assume that the inequality u(— , 0) < u(— , 0) holds in a somewhat uniform way. More precisely, 
for every e > 0, there is a 6 > such that 

u(x, t) < v(y, s) + e if \x — y\ < S and |s| < 6. (8.1) 

This is the case if for example u{— , 0) < v(— ,0) and both u and v are uniformly continuous on 
M. n x {0}. Note that this condition implies that there is a modulus of continuity uj such that 

u £ {x 7 0) < v e (x,0)+uj(e) for all x E R™. 

Where u £ and v £ are the sup and inf convolutions respectively. 

Proposition 8.6. Let u be a subsolution of (|l.ip and v be a super 'solution of (jl.l[) . Assume that 
for every e > 0, there is a S > sucft </ia£ (|8.1|) holds. Then u < v in R" x [0, +oo). 

Proof. Let it e and w e be the sup and inf convolutions, which are a subsolution and a supersolution 
respectively. Since (|8.1|) holds, there is a modulus of continuity u> such that 

u E (x,0) < v £ (x,0) + uj(e) for all x G R". 

We will prove that u £ — v £ < u(e) in K" x [0, +oo). The Proposition follows by passing to the 
limit as e — > 0. 

Assume that supw e — v E — co(e) = a > 0. We cannot assure that the supremum is realized at 
any point, however, for any 6 > 0, there is some point (xo,io) (with to > 0) such that u e (xo,to) — 
v £ (x ,t ) - u>(e) > a - 8. 

Let g be a smooth bump function, for example g(x) = 1/(1 + |^| 2 )- We proceed like in the 
proof of Lemma 13 . 31 we look for the maximum of the function 

w = u e — v £ + 28g(\(x — xq)) 1 

to 

for A <C 1/io- If 8 is small enough, the maximum must be realized at some point (xi,ti) with 
ti > since w(xo,t ) > a + oj(e) but w(x,t) < uj(e) for x large or t large and w < w(e) + 25 on 
t = 0. 

Like in the proof of Lemma 13.21 at the point (xi,ti) where the maximum is realized, the 
functions u e and v £ are C ' (in the sense that they have a tangent paraboloid from above and 
below). So, the equation can be evaluated in the classical sense at that point and we obtain 

d t u e = d t v £ + ^ , 
to 

Vu E = Vi> £ - SWg(\(xi - x )), 
(-A)VV = (-A) 1/2 u e - a(-A) 1 / 2 .g(A(x 1 - x )). 
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From the Lipschitz estimate on u £ and v s , we know that |Vu 6 |, |Vu e < C/e. Let A be the 
Lipschitz constant of H in the ball of radius 2C/e. Then we have 



> d t u e + H(\7u £ ) + (-A) 1/2 u £ 

> 8 t v £ + 5/t + H{Vv e ) + A5X \\Vg\\ LO 

> i - ASX \\Vg\\ L „ - SX {~l\) l/2 g 
c o 



5X 



But this is a contradiction if A was chosen much smaller than 1/to, and we finish the proof. □ 



Existence of solutions 

We now sketch the proof of existence of a viscosity solution to either (|1.2[) or by Perron's 

method. We write the proof for but the same proof would work for (|1.2[) . We choose to 

write it for (jTTTJ) because that case is slightly more difficult since the comparison principle requires 
the inequality at time t — to hold uniformly as in (|8.1|) . 

Let Wo be a uniformly continuous function in R™. We will prove that there exists a continuous 
function u : M n x [0, +oo) that solves and u(— , 0) = uq. 

Perron's method consists in taking the infimum (or relaxed infimum) of the family of all 
supcrsolutions of the equation. There are standard methods, using the comparison principle, to 
prove that this infimum is a continuous viscosity solution. But there is an extra difficulty for 
every particular equation in constructing the appropriate barrier functions in order to prove that 
the infimum of all supcrsolutions is continuous at t = and the initial condition u{— ,0) = uo is 
satisfied. 

Let U be the set of all supersolutions u such that there is some modulus of continuity u so 
that for every x, y £ R™ and t > 0, 

u(y,t)>u (x)-Lj(\y-x\+t). (8.2) 

We start by constructing appropriate barriers to show that the set U is non empty and bounded 
below. 

Let b be a smooth bump function such that b(0) = 1, b < 1 and supp6 = B±. Then, depending 
on the modulus of continuity of Uo, for every e > and xo s K™, there is a 5 > so that 

Uq(x) := b((x - x )/S)u (x Q ) + (1 - b((x - x )/5)) \\u \\ L oo + £> u (x) 
L (x) := b((x - x )/S)uo(xq) - (1 - b((x - x )/S)) \\uo\\ Loa —e< u (x) 

Since Uo andLo are smooth functions, |V/7o|, |(— A) 1//2 f/o|, Vio and |(— A) 1 / 2 io are bounded 
by some constant C . So we can construct a supcrsolution and a subsolution respectively by 

U(x) := U (x) +Ct 
L(x) := L (x) - Ct 

Note that both U and L are uniformly continuous, so we can apply the comparison principle 
of Proposition 18.61 By comparison principle, every u £ IA satisfies u > L for all lower barriers L 
(for all xo S l n and e > 0). Moreover, U is not empty since every upper barrier U belongs to U. 

Let it* be the following function: 

u*(x,t) =liminf inf inf u(y,s) 

r—>Q \x—y\<ru&A 
\t-s\<r 
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For every xq £ 1" and e > 0, L(x,t) < u*(x,t) < U(x,t), since U £lA and for every u £U, 
u > L by the comparison principle (Proposition 18.6]) . Therefore, u* is uniformly continuous on 
R™ x {0} and u,(-,0) = u . 

It can be shown (see [5] for the general method) that u t is a continuous viscosity solution of 
the equation. 
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